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Abstract 



An attempt is made for a new type of analysis of the time-variabihty of the fine-structure 
constant trying to fit the most recent resuh from the laboratory measurements, the Oklo con- 
straint and the data from the QSO absorption lines all in consistency with the accelerating 
universe. 

We have developed a theoretical model of the accelerating universe based on the scalar-tensor 
theory We found that the simplest version of the scalar-tensor theory with an assumed cosmo- 

logical constant A included must be modified in order for the scalar field, serving as the dark energy, to 
provide with its energy density which behaves as a plateau acting like a temporary constant to cause a 
mini-inflation, or an acceleration of the universe. For this purpose we introduced another scalar field x 
providing a temporary potential which traps the original scalar field cr, thus implementing the scenario 
of a decaying cosmological constant; A ~ t^^ in the Einstein conformal frame, realized numerically 
by 10^^^" ^ (10^")^^ in the reduced Planckian unit system, with to ^ 10^" for the present age of the 
universe expressed in units of the Planck time ^ 10^^'^s. This leaves the "fine-tuning problem" solved 
and the "coincidence problem" eased, at least [ll[3l|4]. 

During this trapping process, we naturally expect an oscillatory behavior as a function of 
the cosmic time t. We also showed that this might be observed as a time-dependent fine-structure 
constant. In [5], we tried to fit the most recent results obtained from QSO absorption spectra [B]- 
[2], the constraint from the Oklo phenomenon [TD]-[T3] and the laboratory measurements of atomic 
clocks of Yb+ and IIg+ 15 in terms of the theoretical curves chosen to fit the observed cosmological 
acceleration. 

We obtained a class of fits which we considered to be a reasonable success of the theoretical model 
[5]. We point out that the achieved strong constraint is due to the unique analysis of Single Ion 
Differential a Measurement (SID AM) featuring the measurement at individual redshift [6]- [9], which 
is better suited for comparison with the theoretical assumption of the oscillatory behavior, than other 
approaches making extensive use of the averaging processes |16|-[17|. for which the phenomenological 
analyses tended to produce less constrained results [IS] . 

More important to be noticed currently, however, we have come to know a latest analysis of the 
atomic clocks of A1+ and Hg+ giving the result [TH]; (a/Q;)o — (— 1.6±2.3) x lO^^V^^: which is about 
an order of magnitude smaller than the previous result |15j , and is translated into 




(1) 



where z is redshift and the subscript imphes the present time. 

By extrapolating aja ^ lO^^^y^^ assumed to be time- independent, at this moment, back to the 
epochs ^ lO^^y ago we find Aa/a as small as ^ 10^^. Theoretically, however, today's small value by 
no means implies uniquely the same at earlier epochs [SJ|TB]. In the analysis of [S], on the other hand, 
rather large values of the QSO results can be fitted only by j/q so large compared with ([I]). It 

seems as if we have come to the extent that we may no longer be allowed to identify both of today 
and the Oklo time with the same type of zeros of Aa/a. In spite of certain uncertainties in [6]- [9], 
as well as other results (TBI |T7] , together with still the preliminary aspect admitted in [T5] , it seems 
appropriate and urgent to show if an alternative analysis can be applied basically to the same pursuit 
in terms of the oscillating a and SIDAM, as in [S|. 
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From this point of view, we now assume that we at the present epoch happen to be at the top or 
the bottom of the oscillation, with the Oklo time simply sharing nearly the same height as of today, 
as illustrated in the diagrams (a) and (b) in Fig. 1. Choosing j/q = might be accepted at least as 
a first approximation to deal with the laboratory measurements as small as ~ 10~"'^^y~"'^ [19]. This 
approach, however, requires the separation Atokio ~ 1-95 x lO^y to the Oklo phenomenon to be close 
to one "period" of oscillation. Conversely, the period must be nearly as small as Aiokio. Comparing 
this with the situation in |^ in which the separation between two zeros is not related directly to the 
period, or the half-period, of the oscillation by adjusting the value of j/q, we find ourselves obviously 
more constrained. This poses in fact a non-trivial task as will be also illustrated in (c) and (d) as long 
as we stay in the conventional way of understanding the accelerating universe. 
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Figure 1: Typical plots of the theoretical curve for (Aa/a) x 10^ as a function of redshift z. The Oklo 
phenomenon having occurred « 1.95 x lO^y ago corresponds to z 0.15, indicated in the parenthesis 
of each plot, depending on how the scale factor varies with the cosmic time. The two QSO data 
[6]-[9j are shown; —0.12 ± 1.79 and 5.66 ± 2.67 for z = 1.15 and 1.84, respectively. At the present 
time z = the portion of the "flat" top, or the bottom, is hardly discernible, though it may extend 
over the time-span as long as ^ lO^y, as will be shown toward the end of this article. We commonly 
choose the initial values at ti = 10^" in the reduced Planckian unit system, as in Fig. 5.8 of [T]; 
(Ti = 6.77341501,0-5. — ^tXi ~ 0.21, x']^ = 0, where the prime is for the derivative with respect to 
T = Ini. The lower plots (c) and (d) with /? = arc for the potential in the original theoretical model, 
while the upper plots (a) and (b) are for the Sine-Gordon part of the potential corrected by ([3|) with 
/9 = 0.005, r = 1.2 X 10^"* also with crim as indicated in the diagrams. We applied the technique of 
shifting the evolution variable x as explained in the text on x and 5x. In (c) we cannot make the 
second zero of oscillation to be sufficiently close to the Oklo time as far as f^A and h remain close to 
the observed values J^a = 0.72 ± 0.06 and h = 0.73 ± 0.03 [20l [21], while we manage it in (d) at the 
cost of unacceptably larger values for them, for the reasons we explain in the text. In the plots (a) and 
(b), in contrast, with the potential modified phenomenologically we achieve satisfactory agreements 
with the observational results. 




Figure 2: Computed values of cr, Qa and h are plotted against the evolution coordinate x = log^Q t for 
the modified and the unmodified potentials in the upper and the lower panels, respectively. The same 
initial values as in Fig. 1 were used in obtaining the cosmological solutions. We may re-interpret x 
also as the shifted coordinate x = x + Sx. Today's age of the universe to = 1.37 x 10^°y corresponds 
to xq = 60.204. The arrows marked with the names of plots in Fig. 1 indicate the values xq which 
corresponds to today. fl\ is simply computed in the Einstein conformal frame at xq for x, while h 
is today's value of H in units of i?ioo = 100 km/sec/Mpc = 0.875 x 10"^" in the reduced Planckian 
units, and is estimated hy h = b'{xo)tQ^ /Hiqq, where b' = {d/dT){\na). The values of parameters in 
([3]) are also indicated in the upper panel. 

Let us illustrate how the scalar field a develops with time around the present epoch toward the 
end of mini-inflation. We will also argue how we may re-use the existing solutions of the cosmological 
equation which is highly nonlinear. In the lower panel of Fig. 2, we plot cr as a function of the evolution 
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coordinate x — logj^Q t with the initial value ui = 6.77341501, somewhat different from that used in 
Fig. 5.8 of Ij, but with others remaining the same, imposed at Xi = 10 in the reduced Planckian 
units, so that we find ^Ia and h somewhat smaller than the accepted values 0.72 ±0.06 and 0.73 ±0.03 
pni I21j . respectively, at xq = 60.204 corresponding to today io ~ 1-37 x lO^'^y. We may interpret x 
in Fig. 2 not only to be a real evolution coordinate but also as being "shifted" by a; = a; + 5x. 

In the lower panel of Fig. 2 we pick up a peak at a; = 61.08, for example, chosen for xq corresponding 
to today, as indicated by the arrow marked with (c). We then have 5x — 61.08 — 60.204 = 0.876, also 
allowing us to assign the reasonable set of values JIa ^ 0.77 and h ^ 0.74. The initial time — ti ago 
with ti = 10^'^ when the initial conditions were given originally can be re-interpreted as the occurrence 
at ^0 ~ ^1 — (^0 — ti) + 5t ago, where St/t « dxln 10. This is what we mean by re-using the existing 
solutions. Shifting the initial time may be viewed as another way of revising initial conditions. 

The above choice of xq has been made in such a way that it corresponds precisely to one of the 
tops (or the bottoms) of the curve of cr, hence to the same of Aa/a computed according to 

— = Z-CAtj « 1.84 X lO-^Aa, (2) 
a TT 

which is 4 times larger than in (1) of |5j due to our using more realistic half-spin particles rather than 
spinless particles in the toy model of [U [5], included in the photon self-energy part, but using the 
same regularization method in terms of the continuous spacetime dimension. 

From this "top," chosen to be a zero of Aa for today, we draw a horizontal line to the right 
now in (c) of Fig. 1. The near crossing with the curve of Aa/a is supposed to coincide with the 
Oklo. Obviously we fail because the period of oscillation is too long. We may improve the situation 
by exploiting the observation in Fig. 2 showing that the a oscillation "dwindles" with time in the 
amplitude as well as in the period. Choose another peak at xq — 61.235, for example, further to 
the right in the lower panel of Fig. 2, as shown by the arrow marked with (d). Repeating the same 
procedure as above, we now have a shorter period to cross the Oklo point as in (d) of Fig. 1, but with 
unacceptably larger values of Q\ and ft, as a trade-off. This type of failure appears to be rooted more 
generally, because from a matured "hesitation" behavior the less matter component and the faster 
growing scale factor should follow, as we can easily recognize in Fig. 5.8 of [T]. 

This discussion also suggests a possible way out, faster dwindling of the oscillation as illustrated 
in the upper panel of Fig. 2, perhaps due to the potential "sharpened" toward its bottom. The 
desired potential, partially of the Sine-Gordon type in the original model, is modified according to the 
assumed phenomenological parametrization defined by 

/ \ 2^ 

sm(K(Tj — > sm(K(Tj — p exp 



r 



(3) 



The parameter F in the Gaussian distribution has been chosen to be of the same order of magnitude 
as the range <, 10^"* across which a changes around the present time, as we find in the lower panel of 
Fig. 2. The limiting value (Tim may be defined in principle by Kcrim = (37i'/2) (mod 2tt), but appears 
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slightly different in practice. We consider its local value to be another adjustable parameter. The 
overall factor f3 has been chosen to be 0.005, indicating that the required change of depth of the 
potential is only about 0.5%. With a choice of F = 1.2 x lO^"', (Tim = 45.097392, we pick up the peak 
at xq = 61.0516, as indicated by the arrow in the upper panel of Fig. 2, thus yielding the plot (a) of 
Fig. 1. Note that Q\ and h turn out to fall into the reasonable range caused by a relatively minor 
difference between the behaviors of a in the upper and lower panels of Fig. 2. We point out that 
the fit results in no significant difference in the cosmological behavior of the solution, including the 
equation of state for the dark energy remaining close to —1, basically the same as in Fig. 6 of [S]. 

A similar fit is also shown in the plot (b) of the same Fig. 1. This, though with a somewhat smaller 
h, together with (a) represents our goal in the present article illustrating how successful we can be 
in reconciling the small laboratory value with the large QSO results. It might be worth pointing out 
that it would have been rather easier to find solutions with smaller oscillations, like in the plot (f) of 
Fig. 1 of [5], for example, if the QSO measurements were to yield smaller values. At the same time, 
however, we add that Aa/a could be considerably larger at earlier epochs, like those of CMB or the 
primordial nucleosynthesis, as inferred from Fig. 5.10 of [T] or Figs. 4 and 5 of [S]. 

We are still not sure if the potential modified above only slightly through the phenomenological 
parameters can be derived by revising some of the parameters and initial values of the original theory, 
leaving wider searches for solutions as future tasks. We do not know either why we have i/q ~ 
theoretically at the present epoch. In spite of these questions, we still believe the present attempt to 
be important from wider views. It is also remarkable, on the other hand, to find that the measurements 
at the present epoch have reached the level to affect the assumed earlier temporal behavior of the 
fine-structure constant in a non-trivial manner. 

After all of these discussions comes an obvious message. For further development of the issue, it 
is desperately necessary to start new analyses at least in some of the epochs in which no observation 
or experiment has been available yet with sufficient accuracy. See |22) for an example of insufficient 
rigor of the argument on the constraint expected from the decay -'^^^Re — >-'^*^0s. 

We add a few related comments. We first notice that the correction term ^ applies only to a 
specific cycle of the potential. We may find the similar term which applies to any other cycles for a 
given phase. Let us start with the argument on (Tim mentioned after ([3]), writing 



For the value crim indicated in the plot (a) of Fig. 1 with k = 10, we find 6 = 0.032972. We then write 




(4) 



Ka — k((T — (Tlin) + K(Tl, 



Im — 



v + —{l + 5), 



(5) 



where |w| 1 can be a good approximation. We thus derive 




(6) 
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which can be substituted into ([3]) for an expression to be apphed to any portion of the Sine-Gordon 
potential, as far as we can use a common value of 6. Unfortunately, however, examining the earlier 
mini- inflation to have occurred around x = 26 in the solution used in (a) of Fig. 1, for example, we 
find S = —0.0902, quite different from the above value determined around today. We suspect that a 
may behave differently depending on how it was last trapped by the potential well. In the solution 
under the current study we have only two occurrences of mini- inflation. In all the other cycles, a flies 
so high above the potential that no details near the minima affect the overall behavior of a. 

In Fig. 1 the Oklo constraints were indicated simply by small filled circles, but the real error bars 
are much smaller, nearly invisible. More exact numerical estimates on the fits (a) and (b) in Fig. 
1 result in 0.12 ± 0.01 and 0.006 ± 0.071, respectively, for (Aa/a) x 10^ during the time interval of 
^ 0.1 X 10^ years as an interpretation of Atokio = (1.95 ± 0.05) x lO^y [23]. We refrain from further 
improving the accuracy of the calculation because the present approach does appear to allow its 
implementation, also in view of the current theoretical limitation of the phenomenological conclusion 
from the Oklo phenomenon, as will be sketched below. 

We have the upper bound |Aa/a| <, 10~^ by the conventional Coulomb-only estimate [HIIT^], by 
ignoring the effect of the strong interaction. According to the argument exploiting the near cancellation 
between the two contributions, due to the Coulomb and the strong interactions as elaborated in 
Appendix A of [nil!] however, the process might be interpreted to produce the result which lies, 
roughly speaking, in the range an order of magnitude larger or smaller, representing the size ^ 
0.1 — 0.001 in the scale of Fig. 1, hence leaving the present analysis nearly intact. Model dependencies 
of this type, including asymptotic freedom in QCD hampering the calculation applied to the small 
energy-scale for atomic transitions, are unavoidable before wc find more reliable way to relate the two 
contributions to each other. 

It seems useful to look at a magnified view of the top (or the bottom) portion of Aa/a near the 
present time. An example is shown in Fig. 3 for the plot (a) of Fig. 1. The curve is nearly a symmetric 
parabola assumed centered at xq — 60.204 for computational simplicity. At x — xq = ±1.27 x 10~^, 
we reach the gradient dy/dx — ±7.29 x 10"^, which translates into dy/dz ~ ±3.1 x 10^ corre- 
sponding to both ends of ([T]). The above interval x ~ xq translates also into the time interval 
t w t[j{x a:u)lnlO = ±4.0 x lOV - This is too short to be seen in Fig. 1, yet too long to 

*Tnc second term on the far right-hand side of (A. 3) and on the right-hand side of the second equation of (A. 4) 

ssr^tiJimtfttto- oaiiti[j4ifedcbyam;^oclpaE^el,ei£caii;^rtid[rtfe,ei\grdc^ effectively 

by n^. 
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Figure 3: Magnified view of the top portion of Aa/a = y in the vicinity of today xq = 60.204 
for the plot (a) of Fig. 1. The hmiting values of the gradient ±3.1 x 10~^ in ^ correspond to 
dy/dx = ±7.29 x 10^'' at x ~ xo ^ ±0.000127, or t - to ^ ±4 x lO^y, as indicated by the crosses. 
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